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1/^ , Abstract 

In this paper, we study the nucleon energy spectrum in the ground- 
state (with orbital momentum L = 0) and the first excited state (L = 1). 
The aim of this study is to find the mass and mixing angles of excited 
1^-^ , nucleons using a potential model describing QCD. This potential is of the 

' "Coulombian-I- linear" type and we take into account some relativistic 

effects, namely we use essentially a relativistic kinematic necessary for 
studying light flavors. By this model, we found the proton and A(1232) 
masses respectively equal to (968MeV^, 1168MeV^), and the masses of the 
C^H excited states are between 1564Mel/ and 1607MeV . 



1 Introduction 

In the quark model, the baryons are represented as bound states of three quarks, 
confined by strong forces. Quantum Chromodynamics (QCD) aknowledged as 
the theory of strong interactions, is however not able to describe unambigously 
the strong-coupling regime, which requests alternative theories, like Flux- Tube 
model. Bag model, QCD string model. Lattice QCD, or phenomenological po- 
tential models. We will use in this work a phenomenological potential model 
describing explicitely the QCD characteristics, and taking into account some 
relativistic effects. 

The potential model is essentially motivated by the experiment, and its wave 
functions are used to represent the states of strong interaction and to describe 
the hadrons. The most used is the harmonic oscillator potential which requires 
very simple calculations and qualitatively in agreement with the experimental 
data; on the other hand, the most usual potential models are using non relativis- 
tic kinematics, which is convenient for the heavy flavors systems, but cannot be 
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suitable for the hadrons containing light flavors. 

In this paper, we study the non-strange baryons spectrum in two different 
orbital momentum states within a quark model framework. We calculate the 
mass of baryons in the ground-state (L = 0) with positive parity, and in the 
first excited state {L — 1) with negative parity. The method used is variational 
calculation using Schrodinger equation. The baryon being a three-quarks sys- 
tem, we expose a three-body problem, and the method used to resolve it. 

In section 2, we first introduce the Hamiltonian of our system. We use a 
potential of the "Coulombian-|-linear" form, which reproduces well the confine- 
ment and asymptotic freedom in QCD. For the kinetic part of our Hamiltonian, 
we have added a semi-relativistic correction. 

In section 3, we develop the wave function of our system corresponding to 
ground-state using Jacobi variables. Then we diagonalize the Hamiltonian ma- 
trix and minimise it to extract the baryon mass value. To improve this value, 
we add the hyperfine correction which contains two terms. The first one (called 
the contact term) takes into account the spin-spin interaction, and it correlates 
the splitting of A — iV, S — A in the ground state. The second one, called the 
tensor term, averages to zero for orbital momentum L = Q. 

In section 4, we study the case of P-wave baryons with negative parity. We 
use the same development as for the S-wave baryons to extract masses for the 
first excited state L = 1. In this case, the tensor term is operative, and allows 
us to obtain mixing angles. Our results are then compared with the work of 
Isgur and Karl ppj. 



2 Potential model of QCD 

Several potential models were carried out until now for the determination of 
the baryon mass spectrum in non-perturbative QCD. In general, one uses the 
harmonic oscillator (Q, [3) approach in the quark model . The result obtained 
is in agreement with the experimental data. This approach provides a good 
determination of ground state and excited states of the baryon energy spectrum. 
The Hamiltonian used is in the form: 

3 

Hho = ^ {Ei + VHoinj) + Hhyp) (1) 
i<]=i 

Where is the kinetic part, and the harmonic oscillator potential is in the 
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following form : 

VHo{nj)^lK\V:-7;f (2) 
Hhyp represents the hyperfine correction. 

A more appropriate potential is to be the so-called " Coulombic-I- linear" 
one: complicated non-perturbative effects are assumed to be largely absorbed 
into the constituent quark masses and into a Lorentz scalar linear confinement 
potential, the known short range behaviour of QCD is included in the one gluon 
exchange "Coulombic" potential. In this approach the potential is defined as 
follows P]: 

V{r,,) = -{-^ + -an, + -c)E..F, (3) 

The Fi.Fj factor represents the color term, with a.^ = (Fi.Fj) = — |. 

The constants {as, a, c) are the phenomenological parameters determined 
from an experimental fit. 

In the present work, the calculation of the non-strange baryons spectrum is 
carried out, by using the parameters of reference The results obtained will 
be discussed in sections 3 and 4. 



2.1 Variational Calculation 

The kinetic energy of the relativistic Schrodinger equation describing a system 
with three quarks is in the form : 



E, = ^P^+mj (4) 

But another term, essentially relativistic and more convenient for many-body 
problems than the kinetic part of the Hamiltonian given by equation Q), can 
be used: 

M- in^ 



with the conditions: 



^ 2M,- 2 2AL 



as in reference j7j. Mi represents the quark dynamical masses. They are 
heavier than the quark constituent masses rrii. 



The Hamiltonian of the system is thus written (jH], 



The calculation of the total energy of the system is carried out by the varia- 
tional method. This energy will be minimized w.r.t the variational parameters 
of the test wave function and (Mi). 

The potential V{rij) is treated overall here like a non-perturbative term. In ref- 
erence IH], V{rij) is written in the form of harmonic oscillator potential {Krfj/2) 
plus an unknown U-term (treated like a perturbation) added to shift the energies 
of some states (using the oscillator harmonic wave function) . 

In our work, we concentrate on the determination of the ground-state nu- 
cleon mass, by minimizing the total energy. This method allows us to fix the 
auxiliary dynamical masses (Mi) and the wave function parameters. We use the 
same method to determine the excited-state nucleon mass. The wave function 
parameter and Mj found for L — are different from the ones found ioi L — 1. 



2.2 The three-body problem 

The baryon being a three-body system, so the problem is in the calculation of the 
energy eigenvalue, which contains an integral of nine variables. Before choosing 
the form of the test wave function, one defines the coordinates of Jacobi (7?, A ), 
which make it possible to pass from a three-body system to an equivalent two- 
body system. These two relative coordinates represent respectively the distance 
between the first two quarks ((71. 92), and the distance between the third quark 
qs and the center of mass of the system qiq2- The Jacobi coordinates are defined 
as 12 : 

-> 1 _> _> 

A = — =(ri + r2 - 2r3) 
V6 



Using these coordinates, we obtain the relative Hamiltonian Hn written in 
the following form: 

j/^ = il + il + Af„ + ^ + ^ + ^+y(7?,t) (7) 
lip = Mu {u and d represent the quark flavors) 

P^ 2M„+Md 



With 



The potential V{~p , A ) is written in the form 



A) = --^ 




(8) 



2c 
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3 Application of the model to the S-wave baryon 



In this part, one concentrates on the determination of the ground-states masses 
of non-strange baryons with orbital momentum L=0. The parameter of the wave 
function and masses (M„, M^) are determined using the variational treatment. 

3.1 Determination of the proton wave function 

The total wave function \qqq) of the system is in general constructed from the 
sum of CaYIiX^'^i where (Ca,X)^)^) represent respectively the color, spin, 
spatial and flavor wave functions, with Ca totally antisymmetric. 



The indices A and S mean the antisymmetry and symmetry under the ex- 
change of any pair of quarks of equal masses. 



One restricts oneself here to the spatial wave function for the calculation 
carried out with the Hamiltonian Hr. To calculate the hyperfine correction, 
one will take the total wave function defined in section 4. These corrections 
will give us the mass splitting between nucleon and A(1232). The spatial 
wave function selected to describe our system in the fundamental state 
(^LM — > L = Q, M = 0); is a Gaussian that one develops on the possible states 
of the angular momenta of the systems qiq2 and (73. These states are indexed 
by I Ip, l\) with L = Ip + Ix . The spatial wave function used in our calculation 
is in the following form: 



^lmCp,X) = 5] ^ Q^u-^ipU(W^AmA|LM)xp'''A'^ (10) 



where Ci are coefflcients determined by diagonalisation of the Hn ma- 
trix. The minimization of the system energy allows the determination of the 
variational parameters of the spatial wave function (ap, ax) and the dynamical 
masses (Afu.Md). The average value of the energy over the function ^qq is 
given by the following relation: 



\qqq) = \Color)j^ x \Spatial, Spin, Flavor) 




(9) 



xExp[--{ay + al\^)\Y;^^{np)Y,^^{nx) 



E{ap,ax,Mi) = 



(^ooliJfll^oo) 
(*oo I *oo) 



(11) 
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The only Clebsch-Gordan coefEcients different from zero are those which 
couple the orbital momentum associated to the Jacobian coordinates (7?, A) 
with the total orbital momentum L = 0, and we have: 

\lp-lx\<L = 0<lp + h (12) 

The expression (|12|l yields the constraint {Ip = Ix). So the states of this 
subsystem {~p, A) which contribute to the spatial wave function have all the 
same orbital momentum, starting with Ip = Ix = 0. From parity we have 
P = (— )^ = (— )'p+'^ = +1, so {Ip + Ix) is even. We do our treatment up 
to order 2, i.e. Ip = Ix = 2, since for higher orders the analytical calculation 
becomes very complicated, moreover the contribution of the orbital momentum 
terms ^ > 2 is less important jjl]. The total wave function must obey the 
Pauli Exclusion Principle. In the baryon we have two identical quarks, so the 
spatial wave function must be symmetrical in the exchange of these two quarks 
ri < — > r2, this implies that in the relative system, the function should be even 
in p {~p — -^{^1 — ^))- So in expression (|10|) . Ip must be even. 

Finally, by eliminating the contribution of the odd states in Ip, only two 
terms contribute to the construction of the spatial wave function, which is a 
superposition of \lp — 0,lx = 0) and \lp = 2,lx = 2). So the calculation of energy 
reduces to the calculation of the matrix elements of Hn on the states \lp, Ix) = 
|00) and |22). 

If one notes the spatial wave function in the representation of Dirac, one will 
have: 



|*oo) =ci|00)+C2|22) (13) 

Equation (|13|l shows that the physical state is a mixing of the relative orbital 
momentum states {lp,lx — 0,0) and {lp,lx = 2,2), with mixing coefficients (cl, 
c2). To calculate the energy of the physical state Iv^oo), one starts by evaluating 
the matrix elements of iJ/j on the base {|00), |22)}. The 2x2 matrix is not 
diagonal. If one notes T the kinetic energy and V the potential energy, the Hji 
matrix is written: 



Hr, 



-^1?°°'°°^ \ ^ f (00|(T + V^)|00) (00|V^|22) 



^(22,00) ^(22,22) I (22|F|00) (22 |(r + V^) 1 22) 



(14) 

It should be noted that by analogy with the work of Capstick [5], one has ap 
~ ax , because in the case of the harmonic oscillator ap and ax are proportional 
to the masses m„ and in the case of baryons with the same quark masses. 
So one replaces in the function |5'oo): ctp = c^x = o:. The calculation of the 
elements Ti^i^ in'^ is treated in detail in appendix H6.2I) . The non-zero elements 
are the two diagonal elements Tqcoo and 722,22 • These elements are thus given 
according to (a, Af„,Mrf). 

The calculation of the matrix elements of the potential energy requires to 
evaluate an integral of dimension six on the potential V(7^, A ) which contains 
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Figure 1: Relative coordinates p and A 



With : 



ri2 = V2p (15) 
1 



ri3 = -j=\Jp^ + W6p 



\/3"p . A + 3A2 



r23 = t + 3A2 

and using the form of "Coulombian+ Linear" potential in equation ((SJ 
one calculates now the matrix elements: 

{LM\V{-p,'\)\L'M') ^ Yl {lpmphmx\LM)* {l'pm'/^m'^\L'M') (16) 



X J d-^d\p^''+^'''\^^+^'>^Exp[-a^{p^ + \^)\ V{~p, A) 

xY^^* (17,) y™^* (r!,) (f^,) r^r''- (f^;,) 

To evaluate this integral in the coordinate system (7?, A ) , one 
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uses the hyperspherical coordinates 9) defined as follows 



p = esin(0) (17) 
A = ^cos{e) 

This impHes that : ^'^ = p'^ + X'^ and 9 = arctg{{), < 6* < f 
The integral ^ can be simplified in the following form: 



OO 



With : 

1 



Ki i^is a normalization factor. The problem now amounts to the determina- 
tion of elements Ai^i^ ii'^ and Bi^i^ ii'^. These elements represent the angular 
part in the integral H16(l . The method carried out by calculation of 
and Bi^i^ i'^i'^ is developed in detail in appendix 

3.2 Mass of the proton and mixing coefficients 

The analytic form of the Hamiltonian matrix elements of equation ()14|l is: 



^(00,00) ^ l{_i.7432o+^ + M. + ^ + 2M„ 



+ ^_„44«a + an-L + -|-), (19) 

a Md Mu 

jjmm ^ (0.04396- 0.08154)- (20) 

= (21) 

^(22,22) ^ _o.87160+^+0.5M,+ ^+M„ (22) 

Md Mu 

0.114109 ,1.16667 2.33333a^ , 
+ 0.181177a + — - — + 

a ^ Md Mu ' 



After diagonalisation and minimization of the Hamiltonian matrix, we find the 
following result (using an algorithm of Mathematica 5.0 program) : 

= 1068MeV^ < — {a = 2QQMeV, M„ = Md = 480Mey) 
El^i^ = UUMeV < — (a = 337Mey, M„ = Md = 580Af eV") 
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The physical state energy is the lowest value of these two eigenvalues (-E^j„,-E^j„) 
Table gives the parameters of the potential model considered in this work, 
which were used in The results of the computation of the nucleon and 
A(1232) mass are summarized in table ^ as well as the values of the varia- 
tional parameters (a, M„, Md). These parameters will be used in the calculation 
of the corrections of the spin-spin type in the semi-relativistic model (cf. section 

US. 

It should be noted here that the mass of both nucleon and A(1232) appearing 
in table (0) is the energy of the physical state, which is a mixing of the states 
|00) and |22), where the coefficients of the mixing obtained after diagonalisation, 
are given by: 

|iV, A) = 0.9761 100) - 0.2673 |22) — > En,a = lOGSMeV (23) 

3.3 Hyperfine interaction and calculation of the baryons 
masses 

The potential of the hyperfine interaction is composed of two terms (JSiEl) '■ 

Vhyp = Vc + Vt (24) 

Where the first term (called contact term of Fermi), represents the spin- 
spin interaction between the quarks in the baryon. It is operative only in the 
fundamental state, where the orbital momentum is equal to zero. 

^ o 3 

{ah,Uh) are parameters fitted in reference (see table(l)). The second 
term (called tensor term) represents the static interaction of two intrinsic mag- 
netic dipoles. It is operative only if the orbital momentum is larger than zero. 
It is written in the form: 

N , 

^* = E '^^^WJ^^ {3{S,.r){S,.r) - S,.S,) (26) 

i<j=l ' ^ 

This term enables us to have the mixing coefficients of non-strange baryons 
{S = 0) in the P wave {L = 1). The first term Vc enables us to separate between 
baryons of the same but of different spin. From equations {(O, HSl}, we 
obtain the mass of the nucleon and A(1232) (mass splitting A — N). Table jSJl 
shows the results obtained for the baryons masses, with and without hyperfine 
correction. 
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ah <7h as cr jc{MeV) 


mu{MeV) 


md{MeV) 


0.840 0.700 0.857 0.154 -436 


375 


375 



Table 1: Parameters of the quark potential model [3] 





a 




Md 


Ah 




Exp 




200 


480 


480 


1068 


968 


938 




200 


480 


480 


1068 


1168 


1232 



Table 2: The ground-state energy of baryons (in MeV). Mq is the baryon mass 
without hyperfine correction and Alcor represents the baryon mass with spin- 
spin iVc) correction. 

4 Application of the model to the P-wave baryon 

One carries out now the calculation of the first excited state energy of the baryon 
with L = 1 and negative parity, using the same method as for the S-wave. In 
this case one will build the total wave function, which has to be antisymmetric. 

4.1 Determination of the wave function for L = 1 

In the construction of the spatial wave function, the only possible states for 
(Ip, Ix) which allow to have an orbital momentum L=l and negative parity are 
{Ip = l. l\ = 0) and {Ip = 0, Za = 1) which are noted respectively (^'imj ^im)- 
Note from equation © that 'i'ljij is even under the exchange of the first two 
quarks, while the analogous wave function 5'^^,/ is odd, since p and A respectively 
indicate mixed antisymmetry and mixed symmetry under this transposition. 
These states form a representation of dimension two of the permutation group 
5*3, which allows to exchange each pair of quarks in the baryon. 
These two functions have the following form: 

Cio ^1" i^mpOmx \ IM) 
xpExp[-^a^p' + X')]Yr^{np)Y,"^^{nx) (27) 

Coi J2 ^01 {Omplmx \ IM) 

mp,mx 

x\Exp[~\a\p^ + \^)]Y^^{np)Y^^{nx) 

One will have also to consider the coupling of the spins S with the orbital 
momentum L to build the total angular momentum J = L + S . In short, here 
is the construction of J: 



*^m(7?, A) = 
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The spin of the quark system : S* = i, | , is constructed as the triple product 

One has L — 1, one thus distinguishes two constructions of J for the doublet 
and the quadruplet of spin : 



'2)S=lL 



1 



2 ' 2 ' 2 



The nucleon states of total angular momentum ^ = ^ , § are a mixing of the 
doublet and quadruplet of spin, for example: 



) = Oil 




+ /3i 


) = ^2 




+ /32 



(28) 



Thus, the wave function will be built as 



'Pi 

2 



(29) 



The construction of the spin-flavor wave function is carried out in analogy with 
the notation of Karl and Isgur and collaborators (0, |Hj)- The flavor-mixed 
antisymmetry and symmetry combinations of "uud" are: 



p 



—= {udu — duu) 
v2 



(30) 



— — —j= {duu + udu — 2uud) 



In the same way, the spin states are: 



x| = ^(inT)-UTT)) 

xj = -^(UTT) + ITiT)-2|TTi)) 



(31) 



One must add the spin wave function for the quadruplet which is completely 
symmetric: 



xf = ITTT) 



(32) 



Thus, the total wave function will be built according to the following com- 
binations: 
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Z Z / Z 2 2 2 2 



A^'^Af(2 '2 7^ CAxl^i'^'N'^'lM + '^N^lM) (33) 

A^^M(i^^")) = C^'i>i^(*?Mx| + *W|) 

From equation H33() . one sees that we have two groups of states: N* and A*. 
Each group has its own flavor wave function. 

4.2 Excited nucleon mass 

For the P-wave excited-state nucleon, the Hamiltonian matrix is written as 
follows: 

/ m Hn \Q1) {01\Hnm \ 

{lQ\Hn\Ql) (10|Hfl|10> J ^-^^^ 

By the same method used in calculation of the baryons masses in the ground- 
state L = 0, one finds the elements of the Hb matrix for L = 1: 



1 7372 
(01|Hfl|01> = -(-1.7432-1- 2.7A45a+-^ + Md 

2 a Aid 

2/1 2 , 2ml. 
+a ( \ ) H -) 

Hid mJ M^' 

(01|i/fl|10) = {lO\HB.m 

= (0.0439 -O.OSlSa^)- 

a 

0.1141 0.1812 1.1667a2 Md 
\ h — - + — 

a a Md 2Md 2 

2.3333^2 ml ,^ 

After diagonalisation of this matrix, one obtains two eigenvalues of the 
Hamiltonian Hb ■ 



(10|fffl|10) = -0.8716 
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1 1, mf, 2mi 0.2027 
H], = - -1.7432 + — ^ + 2M„ + 1.0745a 

,1.3333 2.6666, , r „ „o.o 0.0310 ,0.0193 0.0193, 
H^T- + " - {-0.2863 + — - — — + -r^ a 



Ir . to3 , , 2m2 0.2027 



Hi = -[-1.7432+ —^i + ^ + 2M„ + ^^^^^ - 1.0745a (36) 
2 Md M„ a 

,1.3333 2.6666, , r „ „o.o 0.0310 ,0.0193 0.0193, 
+ (^r^ + -Tl—)^ + {-0.2863 + — - — — + -r^ a 

Md Mu a- Md Mu 



-0.7850a2 - (- + )a^ + ( ' ^ + ' , + )a^}3 

^ Md Mu ' ^ MJ Ml MuMd' ^ 



The minimization of these two eigenvalues yields the following result 



Hi 1590 MeV 

Hi 1755 MeV (37) 

The minimal value (1594 MeV) corresponds to the nucleon mass in the first 
excited state L = 1 for the values (a = 192MeV,Mu = Md = 1485Mey). In 
order to separate between the nucleon states of different total angular momen- 
tum, one will apply the hyperfine potential defined in equation H24|) which is 
composed of two parts, the contact term and the tensor term. 



4.3 Hyperfine correction 

The orbital excitation i = 1 of the baryon is carried by only one pair of quarks, 
while the two other pairs have zero orbital momentum. This decomposition 
of orbital momentum enables us to note that the tensor part of the hyperfine 
interaction is operational only in the quark-pair with L=l. On the other hand, 
the two other pairs with orbital momentum equal to zero are controlled by the 
contact term of the hyperfine interaction. One will calculate first the contact 
part, then the tensor part. 



4.3.1 The elements of the contact matrix 

The following calculation was carried out for only one pair of quarks, then the 
result was multiplied by three. Furthermore the result does not depend on the 
total angular momentum. Recall that the total spin for our system takes the 
value (5 = ^, 5 = |) and thus we are left to calculate the elements of the 
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following matrix : 

/ {N^Pm{J)\V,\N^Pm{J)) {N^Pm{J)\V,\N^Pm{J)) \ 

(7V2p^,(j)| V, \N^Pm{J)) {N^Pm{J)\ Vc \N^Pm{J)) 

\ 

Using equation one obtains: 



(38) 



(N^Pm{J)\Vc\^^Pm{J)) 



Z{N^Pm{J)\V}''\N^Pm[J)) 
3 



3 %i^ah 2 a\ 
2 3M2 3 
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S'i.S'2 



(TV^FmI J)| K \N^Pm{J)) = 3 (iV'FAf( J)| K'' \N^Pm{J)) 

_ 3 87ra,, 2 al 



2 3M2 3 



((' 



S1.S2 



4) + (4 



s'i.s'2 

A 



x{{^>lM\Exp{-alrl^)\^>l,t) 
+ {^%,\Exp{-<jlrl^)\^'[,,)) 



(39) 



(40) 



After integration of the spatial part, and using the Jacobi variable change 
(ri2 — V2/o) one obtains the following result : 



(41) 



{^',M\Exp{-<jlrl,)\^P,,) 



{'f'^,,\Exp{~2alp')\^P,,,) 



(a2 + 2a2)f 

Equations { (p!^ . (001)} are written then in the following form: 



(42) 
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{N''Pm{J)\V}''\N^Pm{J)) 
{N^Pm{J)\V}^\N^Pm{J)) 



(43) 



{N^Pm{J)\V^^\N^Pm{J)) 



a 



where D — ^§jk- The result of equation H43I) docs not depend on the total 
angular momentum J. The non-diagonal elements of the V^'^ matrix are equal 
to zero. For the case of the A*, the tensor part of the hyperfine potential does 
not contribute to the correction of the mass, however the contact term gives: 



{A^Pm(J)\V,\A^Pm{J)) 



(A2Pm(J)|K''|A'Pm(J)) 
3 



{ 



3 2ah Sal 



2 3Af 2 30F 




3 8Dal 



2 3^ \ 4(^2 + 3^2)1 4(a2 + 3a2)f 



3 SDa'i 1 



2 3V^ 4(a2 + 3a2)i 



-3 



(a2 



(44) 



4.3.2 The elements of the tensor matrix 

The interaction between magnetic dipoles appears in the case L = 1. The 
elements of the tensor matrix depend on the total angular momentum J. Using 
the formulas of reference |12j . we find the following results : 



N^Ps 



N^P. 



N^Ps 



Vt 



12 



12 



iV2p3 



N^Ps 



3D 
~2 



1) Ml) 



3/5 
V2 

= 



N^Ps 



D 
"2" 



(45) 



15 



12 



2 
2 



2 

15 D 
T 2" 





12 



N'^Pi 



{n,\p-\icosH,-l)\^{,) 



with : 



(*fi|p-3(3cosX-l)l*?i) = 



q;2 

15 



(46) 



4.3.3 Energy spectrum £ind mixing angles of the first excited states 
of the baryons 

After having calculated the matrix elements of the hyperfine correction, the to- 
tal matrix is then written as follows: 



For J = i: 



yhyp — D 



a 

(a2+2<T^) 



For J = |: 



Vhyp = D 



1 

/lO 



(47) 

\ 

/ 

(48) 

After diagonalisation of these two matrices, one finds the following result (in 
MeV): 

For J=\: 



1 

\/lO 



(a2+2<T^) 



(ifij + RhW) 



2 2 



1590 - 26 
1590+ 10 



2 2 



(49) 



For J= |: 
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For the A* one obtains : 

(Hr + Hhyp)\A>= (1590+ 17)\A> (51) 

A physical state of total angular momentum J is defined as a mixing of states 
of spin momenta {S = ^ and |). The mixing angles are defined in the following 
equation, see ^31 ■ 

\N}{min) >=-sm9 \^Pj > +0080 \^Pj > 

\N}{max) >= cose \^Pj > +sin6' fPj > (52) 
This gives the following result : 

iNUmin) 0.512 I^Pi > +0.859 PPi > 

2 2 2 

|7V|(TOaa;) >= 0.859 |*Pi > -0.512 pPi > (53) 



For J = i 



From these equations we obtain the mixing angle noted 9s 



tan(9, = -^^^ =^es = -30.8° (54) 
0.859 ^ ' 



For J = I : 



iNtimin) -0.107 I'^Ps > +0.994 pPa > 

\Nt{max) >= 0.994 j^Pa > +0.107 pP| > (55) 
From these equations we obtain the mixing angle noted dd '■ 

tanOd^^^^ =4> 6*^ = 6.1° (56) 
0.994 ^ ' 



Without hyperfine correction we have found the mass of the excited nu- 
cleon {L = 1) equal to {M^ — 1590Mey) for the following values {a = 
192MeV, Mu ^ Md 1485Mey). 
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M' 

cor 

"Coulombian + linear" 


M' 


ax 

"Coulombian + linear" 


ax 
1 


7Vt(1535) 


1564 


1490 


0.859 


0.850 


7V1(1650) 


1600 


1655 


-0.512 


-0.530 


7V|(1520) 


1573 


1535 


0.994 


0.990 


7V|(1700) 


1620 


1745 


0.107 


0.110 


At (1620) 


1607 


1685 


1 


1 


At (1700) 


1607 


1685 


1 


1 



Table 3: Masses M'^^^iva MeV) and mixing angles Qx of baryons for L=l with 
hyperfine correction calculated by "Coulombian + linear" potential model and 
compared with Isgur and Karl result 0. 



Table (3) contains the nucleon masses {M'^^^.) and the mixing angles (ax) 
with hyperfine correction for "Coulombian + linear" potential model, compared 
with the work of Isgur and Karl U for L = 1. 

The hyperfine correction to the energy spectrum of the excited nucleon is small 
compared to that found in j^. 

We remark that for the nucleon with total angular momentum J = i , one finds 
from equation H54|l a mixing angle 9^ — —30.8°, which is very close to the result 
of ^ 6*5 = —31.7°. For the nucleon with total angular momentum J = |, one 
finds from equation H56|l a mixing angle 9d = 6.1°, which is very close to the 
result of n 6*^ = 6.3°. 

Equations H47I I48|l allow us to have the mixing angles ax ■ They depend on the 
spatial wave function parameter a and the potential parameter ah contrary to 
the work of Isgur and Karl |12) . where the ax are independent of any choice of 
parameters. 

5 Conclusion 

The choice of the potential used (Coulombian-|-linear) in our model and the ad- 
dition of the semi-relativistic correction of the kinetic part of the Hamiltonian 
allowed us to obtain a good result for the masses of the non-strange baryons 
in S-wave. We separated the nucleon state from the A(1232) using the contact 
potential {Ma — Mn — 200MeT^). This difference in mass between the nu- 
cleon and the A(1232) is small compared to that found by Isgur and Karl J3] 
(Ma - Mn = 260Mey). 

For L = 1 we obtained the excited nucleon mass without hyperfine correc- 
tion (Mjv* = 1590Mey). The result is very close to that of reference 
On the other hand the hyperfine correction enables us to separate between 
the states of the excited nucleon. This correction is very small compared 



18 



to the result of pp. The difference is due to the dynamical mass used in 
the hyperfine correction which is heavier than the quark constituent masses 
(M„ = Md = ASOMeV > to„ ^ ma = il^MeV for L = 0) and (M„ = Md = 
1485MeT^ > m„ = = 375Mel/ for L = 1). 

For the case of the mixing angles we obtained the same result as found by Isgur 
and Karl PP, on the other hand our elements of the hyperfine potential matrix 
depend on the parameter a of the spatial wave function and the parameter a of 
the linear part of the potential. 
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6 Appendix 

6.1 Wave function form and development 

One takes a function of the "Gaussian development" : 

*LM(7,t) = Y,Y.Ci^i^.Ni^i^{lpmplxmx\LM)xp^^\^^ (57) 
X Exp[- \ {ay + alX' Wl^" (Jl,)y-^ (Q,) 

{IpTriplxnix I LM) arc Clcbsch-Gordan coefficients (one couples the orbital 
momenta associated to the Jacobian variables (7?, A ) with a total momentum 
t = l^+h )• 

In the case of baryons with at least two identical quarks, the above expression 
of I'^lm) must satisfy the constraints imposed by the Pauli principle. If the 
two identical quarks are numbered 1 and 2, their distance is expressed by the 
Jacobi variable "/? , and the spatial wave function must be even in p. Moreover 
the parity of the proton is positive, and according to the following relation 
(P = (—1)^), one deduces that the total orbital momentum must be even. 
Parameters of the spatial wave function {ap,ax) are given by minimizing the 
total energy E of our system (^). For that, one will calculate the kinetic 
energy and the average potential energy of the system. 

6.2 Calculation of the average kinetic energy 

The average value of our kinetic energy is written as follows: 



E E ^i.i>^^i'A {lpmplxmx\LM)* (58) 



{l'pm'/xm'x\L'M') x j crpdXp^''\^ 



xExp 



--{ay + aiX^) 



1^7''*(n,)F™^*(0,) 



p2 p2 

2/Up 2/iA 



P^'pX^'^ 



X Exp 



1/22, 2 \2\ 



With 
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lp + l^+3 

'"'^ ~ Vr(i+;p)r(f+;;,) 



P' 



1 J. /y,2^\ 

dr \ dr J 

1\ T' — 1\ T ^ ^ Mi 4- ™d 

and : 



- 2 2 

-a r 



2Mu+Afd 

(;^(?^ + 1) - a^{2lr + 3)r2 + aV^) 
1 



xr ''Exp 



2 2 

-a r 



We can now resolve equation H58|l as follows : 



With: 



(59) 



M" - L^] 



RiPn 



(60) 



R (Pi) = ^ / ^^^^ (^'"+'^+'£^2;^ [-a>2] 



xX^^+^'^Exp [-alX^] Sl^l'^Srn^m'^Si^yJ„i„n', 



(61) 
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Exp[~ay]) (63) 



X A 



6.3 Calculation of the average potential energy 

As a first approximation, one takes (ctp — ax = a) . The "Coulombian + linear" 

potential is written according to the variables of Jacobi in equation (jS]) . 

In order to calculate the average potential energy Vi^i^j' i'^, we are going to 

determine the following coefficients [Ai v ,Bi i^ji v ) which represent the 

angular parts of this integral. Using the hypersphcrical coordinates {^^ ,6^, 
one finds the following result : 



^p^X 



1 


V3~X +-p 


1 


V3~X -~p 










1 


- H 


1 


1 

V2 




1 






'in,) 







(64) 



(65) 



Where : 

dV,, = dQpdflx sin^ flcos^ 9d9 , dfii = sm9id9id(pi 
KC' ("5) = A,;, sin'- 9 cos'^ y"" {%) y™^ (f7,) 



For that we use the following definition : 



2 

i= 1 d^sin^cos^ (0) (0) + i (/+ + /-) 



(66) 



V2 y dflsin^cos^ 9^1,^ (9) if,^,^ (9) + ^ (/+ + la) 




(67) 
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With: 



= ldn,\H±H\-'Y,X{^5)Y,^il{^5) (68) 

"2 

= JdO sin^ e cos^ Oipt^,^ {6) ip,,^,'^ (0) J dnpY;^'* (O^) Y^T'' (Op) 



l—O m— — l 

= J dn5\n±i^\Y^'*{n5)Y^' {n^) (69) 

d9 sin^ cos2 0ipl,^ (0) ip,,^,^ (0) I dilpY;^'* (Op) F,"^ (Op) 



X /do,y™-(o,)r<(o,)f; £ (^)'a,-i^rr(^^p)i^r(^^A) 



Where : 



ifii.i^ (0) = Di,i^ sin'" sin'^ 61 (70) 



D : is the normaUzation factor. 



"'^'M^i^-^rri^ ^^^^ 

( r< = niin(ri,r2) , r> = inax(ri,r2) ) . 
From our choice of potential ( Coulombian + Linear ) we have : 

Expanding the two parts and 7^, we obtain : 

uiin[lp+l'^,lx+l'x] 

Ib= E (T)' (-)'-"M-r^ Ml [61 + 62] (73) 

/=0 

It= E (T)'(-)'""M-r^MiM2X [ai+a2] (74) 
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With 







(75) 



Ma = /^+'-'^ ^{2lx + 1) (2Z^ + 1) 







(76) 



arctan(2) 



sin' 6* 



COS' 



I+IQ 



(77) 



62 = J de sin^ cos^ e^Y^* (6) (e) X (2 X 3') ^ 



1 / cos' 9 



sin'+i 6* 



arctan(2) 



(78) 



arctan(2) 



ai = 



/ 



rf^sin^ ^cos' 0^Y"* (^) *L''^ (^) 



(79) 



sm 



i+2 , 



2x3'+V (2Z + 3) cos'+i 6* V2 x 3'-V (2Z-l)cos'-i 



sm 



012 



2 

/ 



'sin^ ^cos^ 0^^£^'* (9) ^^l!''* (6) 



(80) 



arctan(2) 

'3i+2\ f 



COS 



,i+2^ 



{21 + 3) sin'+i 61 



3'y 2'-icos'^ 
2"; (2Z- l)sin'-i(9 



Finally the total average potential energy is the sum of matrix elements 
Vi^i^^i'^i'^ multiplied by Clebsch-Gordan coefficients : 

Vlm,l'M'= Y1 E {lprnplxmx\LM){l'^m'/xm'x\L'M')xVi^i,,i,^r^ (81) 



With : 



L = Ip + Ix , 



(M = mp + toa) 
(M'=m;, + m'J 
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The matrix elements Vi^i^^i^i'^ lead to enormous calculus, so we are in the 
obligation to restrict ourselves to the lower orbital excitations (pratically we 
limit ourselves to {IpJxJ'pJ'x < 2) ). Doing such restriction leads to a very 
simple analytically expression wich can be solved. 



6.4 Some relations 



I 



(82) 
(83) 



diiY/^'* {n)Yi"' {n)Yii 



(84) 



(2;i + l)(2Z + l)(2;2 + l) 



^1 



An 
I h 



-mi 7712 



h I h 





E 



h k L 

7771 ^2 — Af 



h h L' 

7771 "l2 — M' 



1 



2i+ 1 



JLL'OmAI' 



(85) 
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